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Differential elliptic flow and particle spectra are calculated from covariant Boltz- 
mann transport theory taking into account the finite transport opacity of the 
gluon plasma produced in Au+Au at E cm ~ 130 A GeV at RHIC. The solutions 
are shown to depend mainly on the transport opacity, \ = J dzatPg- The elliptic 
flow saturation pattern reported by STAR indicates that Xb—o ~ 25, i.e., the gluon 
plasma is ~ 80 times more opaque than the pQCD estimate based on HIJING. 
Such large opacities are also consistent with the measured charged hadron spectra. 



1 Introduction 

At present we have a very limited understanding of the properties of the par- 
tonic environment created in heavy-ion collisions at RHIC. Predictions for the 
density of the gluons produced vary by a factor of five depending on the model 
considered. Elliptic flow, V2(p±) — (cos(20)) px , the differential second mo- 
ment of the azimuthal momentum distribution, and the high-p_L s^pression 
of the particle spectra have been the subject of increasing interestElu'Erua be- 
cause they provide important constraints on the density and effective energy 
loss of partons. 

The origin of the remarkable saturation of elliptic flow iWpjJ — * 0.2 
above p± ~ 2 GeV reported by STAR at Quark Matter 2001B is an open 
question for theory. Calculations based on inelastic parton energy lossa do 
predict saturation or decreasing V2 at high p±. These calculations are valid 
for high p±, where collective transverse flow from lower-pj^ partons can be 
neglected and Eikonal dynamics is applicable. However, a constant spatial 
anisotropy was assumed throughout the evolution, while in reality, it decreases 
and probably even changes sign. This is likely to reduce the generated elliptic 
flow much belowu the preliminary data. 

Ideal hydrodynamicsEl, the simplest theoretical framework to study ellip- 
tic flow, agrees remarkably well with the measured elliptic flow dataH up to 
transverse momenta ~ 1.5 GeV/c. However, it fails to saturate at high p± > 2 
GeV as does the preliminary data. 

A theoretical problem with ideal hydrodynamics is that it assumes local 
equilibrium throughout the whole evolution. This idealization is marginal for 
conditions encountered in heavy ion collisions!!!. Covariant Boltzmann trans- 



port theory provides a convenient framework for nonequilibrium dynamics 
that depends on the local mean_fcee path X(x) = l/an(x). 

Parton cascade simulationalu show on the other hand, that the initial 
parton density based on HIJINGQ is too low to produce the observed elliptic 
flow unless the pQCD cross sections ace artificially enhanced by a factor ~ 
2 — 3. However, gluon saturation modelsa predict up to five times higher initial 
densities, and these may be dense enough to generate the observed collective 
flow even with pQCD elastic cross sections. 

In this study, we explore the dependence of elliptic flow and the high- 
p± suppression of the particle spectra on the initial density and the elastic 
gg cross section. Though parton cascades lack at present covariant inelastic 
energy loss, elastic energy loss alone may account for the observed high-p^ 
azimuthal flow pattern as long as the number of elastic collisions is large 
enoughu. 



2 Covariant parton transport theory 

We consider here, as in the simplest nonlinear form of Lorentz- 

covariant Boltzmann transport theory in which the on-shell phase space den- 
sity, evolves with an_elastic 2^2 rate. We solve the transport equation via 
the MPC algorithrn3 r Jjdaich maintains Lorcntz covariance using the parton 
subdivision techniqueE2rE3. See Ref. u and references therein for details. 

For a given nuclear geometry and formation time, the solutions of the 
nonlinear transport equation has been shownu to depend mainly on the trans- 
port opacity X = J dza t Pg = -/V(sin 2 9 cm ) and the impact parameter b. Here 
o~t is the elastic transport cross section, N is the average number of collisions 
per parton during the whole evolution, while 9 cm is the collision deflection 
angle in the cm. frame. To good accuracy the transport opacity factorizeaj 
as X = C(b)a t {T )dN g /dn. 

We label our results by the transport opacity x an d impact parameter b. 
Furthermore, we quote our transport opacities relative to that for the pQCD 
minijet gluons predicted by HIJING [dN/dr] = 210, o gg ^ gg « 3 mb => er f ss I 
mb, Xb=o « 0.3]. 



3 Numerical results 

The initial condition was a longitudinally boost invariant thermal Bjorken 
tube at proper time to = 0.1 fm/c with uniform pseudo-rapidity distribution 
between \r)\ < 5 and transverse density distribution proportional to the binary 
collision distribution for the two gold nuclei. Based on HIJING, the pQCD 
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Figure 1. Impact-parameter-averaged elliptic flow as a function of p± for Au+Au at 
y/s = 130A GeV for different transport opacities with hadronization via local parton-hadron 
duality (left) and independent fragmentation (right). 




Figure 2. Charged hadron elliptic flow vs. centrality (left) and negative hadron pj_ spectra 
for b = 6 fm (middle and right) as a function of the transport opacity for Au+Au at 
y/s = 130A GeV. Left and middle plots are for hadronization via parton-hadron duality, 
while right plot is for independent fragmentation. The spectra are normalized to 1 at = 2 
GeV. 

jet cross section was normalized to yield dN g /drj = 210 in central collisions 
and the initial temperature was chosen to be Tq — 700 MeV. 

Two different hadronization schemes were applied. One is based on local 
parton-hadron duality, where each gluon is assumed to convert to a piom. 
The other hadronization prescription is independent fragmentation, where we 
considered only the g — * tt^ channel. See Ref. 13 for details. 

Fig. [l] shows the impact-parameter-averaged elliptic flow as a function 
of p±. With increasing p±, elliptic flow increases until p± ~ 1.5 — 2 GeV, 
where it saturates. To reproduce the preliminary STAR data, ~ 80 times 
more opaque initial gluon plasma is needed than the pQCD prediction from 
HIJING. Surprisingly, the results show no sensitivity to the applied hadroniza- 
tion prescription. 

Fig. 2 shows that such large transport opacities are also consistent with 



the preliminary charged hadron spectra measured by STARllJ. Hadronization 
via parton-hadron duality yields too little suppression at high p± because it 
only incorporates quenching due to elastic energy loss. However, with the 
additional quenching due to independent fragmentation, the parton cascade 
results approach the preliminary STAR data. 
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